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We theoretically investigate the competition between charge-ordered state and Mott insulating state 
at finite temperatures in quarter-filled quasi-one-dimensional electron systems, by studying dimerized 
extended Hubbard chains with interchain Coulomb interactions. In order to take into account one- 
dimensional fluctuations properly, we apply the bosonization method to an effective model obtained by the 
interchain mean-field approximation. The results show that lattice dimerization, especially in the critical 
region, and frustration in the interchain Coulomb interactions reduce the charge-ordering phase transition 
temperature and enlarge the dimer-Mott insulating phase. We also derive a general formula of the Knight 
shift in the charge-ordered phase and its implication to experiments is discussed. 

KEYWORDS: charge-ordering, quarter filling, quasi one-dimension, extended Hubbard model, lattice 
dimerization, Mott insulator, molecular conductors 



Since the experimental observations of the charge-ordered 
(CO) state in quasi-one-dimensional (ID) organic 2:1 com- 
pounds, 1 e.g., (DI-DCNQI) 2 Ag 2 and (TMTTF) 2 X 34 (X: 
a monovalent counter anion), in which the charge distri- 
bution becomes disproportionated in a regular way due to 
the electron-electron Coulomb repulsion, these systems have 
been extensively studied, experimentally 5 as well as theo- 
retically. 6 In (TMTTF^AT, the CO state is easily destroyed 
by external pressure 7 - 8 but nevertheless the system remains 
insulating, which is due to a competing state, i.e., the so- 
called dimer-Mott insulating state where the charge local- 
izes on each lattice-dimerized bond. On the other hand, in 
(DI-DCNQI)2Ag, a recent X-ray structural analysis 9 revealed 
the existence of a three-dimensional mixtured pattern of CO 
and bond dimerization in the "CO phase". In this paper, we 
theoretically investigate such interplay between the CO state 
and the lattice dimerization based on a quasi-one-dimensional 
model, especially their properties at finite-temperatures. 

Previous theoretical studies on the CO phenomena in these 
quasi- ID materials have mainly been devoted to ground-state 
properties of the ID quarter- filled extended Hubbard model 
(EHM) with the on-site and nearest-neighbor Coulomb repul- 
sions U and V as a minimal model. 6 In this model, CO in- 
sulator (COI) is stabilized at T — in the large U and V 
region, while the system is in the metallic state described by 
the Tomonaga-Luttinger liquid (TLL) otherwise (for positive 
U and V). 10-12 When the lattice dimerization is included by 
alternation in the transfer integrals, the TLL phase is trans- 
formed to the dimer-Mott insulating phase, and this competes 
with the CO state. 1315 There, the fluctuation effects in the 
competing region have been focused based on the bosoniza- 
tion theory 16 and critical fluctuation appears on the boundary 
between the two competing states 17-19 where its universality 
class is the two-dimensional Ising-type. 

However, these results cannot be directly applied to the 
analysis of the actual materials, since purely ID electronic 
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Fig. 1 . Schematic illustration of the model investigated in the present study. 
Each ellipse denotes a molecular site and the a (b) axis shows the chain 
(interchain) direction. The thick (thin) lines along the a-axis shows the 
lattice dimerization of the stronger (weaker) bonds. The solid and dashed 
lines connecting the chains represent the interchain interactions, V± and 
Vj_ , respectively. 



models do not show any phase transition at finite tem- 
peratures. In our previous work, 20 in addition to the ID 
quarter-filled EHM, we have taken into account the interchain 
Coulomb interaction V± and investigated the CO phase tran- 
sition whose transition temperature Trjo becomes always fi- 
nite due to the dimensionality effect. We found that the inter- 
chain interaction can transform the system from TLL metallic 
at V± =0 to CO insulating by infinitesimal V± ^ in the 
critical region and thus the interchain interaction greatly en- 
hances the COI phase in the U -V plane. 

In the present study, we consider a system shown schemat- 
ically in Fig. 1, which is adopted from the crystal structure 
of (TMTTF^AT. This is an extension of the model studied in 
our previous work 20 to include the lattice dimerization along 
the chain direction, as well as different Coulomb interactions 
in the interchain direction. As in Ref. 20, in order to take 
full account of the ID fluctuation effects, we treat an effec- 
tive ID model by the bosonization theory, by applying mean- 
field approximation to the interchain interactions. Our main 
results are summarized as follows, (i) By deriving an analyt- 
ical expression of Tqo an d computing it by the renormaliza- 
tion group (RG) approach, we find that the lattice dimeriza- 
tion largely suppresses the COI phase and results in the dimer- 
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Mott insulating phase, especially in the critical region, (ii) The 
different interchain Coulomb interactions act as geometrical 
frustration for the CO state, therefore the COI phase is sup- 
pressed when these are comparable to each other, consistent 
with previous works. 21 (iii) A general formula for the Knight 
shift below the CO phase transition temperature is derived as 
a function of the magnitude of order parameter. 

We start with the model Hamiltonian -/Jqid = + 
H±, where H 3 ir> represents the jth extended Hubbard chain 
with finite lattice dimerization 6 4: 



+ U ^ n itj> i n itj>l + ^ E n itj n i+1<j 



(1) 



where c\ ■ s is the creation operator of an electron with spin 
s(=T/l) at the ith site on the jth chain. The density operators 
are n iJ)S = ct J s c iJiS - \ and mj = n iJA + n iJd . The 
second term H± expresses the interchain interactions, given 
by 



riij' + V{_ E n h3 n i+i,]'i 



(2) 



where the V^-term denotes the interchain interaction between 
electrons on the nearest-neighbor sites (along the solid hori- 
zontal lines in Fig. 1) and the V[-term denotes the interaction 
for the next-nearest-neighbor sites (along the dotted diagonal 
lines). When V± = Vj_ = the model reproduces the above- 
mentioned ID dimerized EHM. 

As in Ref. 20, we assume the Wigner crystal-type CO pat- 
tern with a 2-fold periodicity along the chain direction, 22 
and apply the interchain mean-field approach 23 to H± [eq. 
(2)]. The two interchain Coulomb interactions favor differ- 
ent CO along the interchain direction, namely, V± favors the 
antiphase pattern, (nij) = (— l) r+ ^n, while V[_ favors the 
inphase CO with (riij) = (— l)*n, where n(> 0) is the am- 
plitude of CO to be determined self-consistently. This results 
in geometrical frustration arising from a competition between 
two types of CO states. 21 The resultant effective ID Hamilto- 
nian can be written for both cases as 

H = 1 E t 1 + (-i)'^] (4,sC i+ l,s + fc.c.) 

i,s 

+ u E " 4 >t + y E ji ' Ui+i 



(3) 



where the chain index j is omitted and N is the total num- 
ber of sites in a chain. The coordination number of adjacent 
chains, which is given by 2 in the system of Fig. 1, is writ- 
ten as z. The effective interchain interaction V± is given by 
V± = V± — V^_, then the interchain frustration simply reduces 
the effective interchain interaction in the present scheme, 
therefore acts as destabilizing the CO state, similarly to pre- 
vious studies on different models. 21 

In order to take fully into account the ID fluctuation ef- 
fects, we treat the Hamiltonian (3) based on the bosoniza- 
tion theory. First, we focus on low-energy excitations near 
the Fermi wavenumber ±kp — ±7r/(4a) (a: lattice spacing), 



by integrating out the high-energy states. 11 Next, following 
the standard bosonization theory, we express the Hamiltonian 
in the bosonic phase variables and separate it into the charge 
and spin parts; H = H p + H a . The spin part H„ takes the 
same form of the effective theory of the ID isotropic S = 1/2 
Heisenberg chain, in which the coupling constants depend on 
(5d as well as n. This indicates that the spin excitation is gap- 
less in the both insulating states, and the competition between 
the COI phase and the dimer-Mott state is characterized solely 
by the charge part of the Hamiltonian, H p . 

Let us first focus on the CO transition temperature Tqo- 
Although the phase Hamiltonian is derived for the full order 
of n by the above operation, in order to evaluate Tqo the ex- 
pression up to the first order in n is sufficient, since a recent 
numerical study 24 shows that the transition is continuous. In 
this order, the charge part is given by H p = J dx(H p o + H p i ) 
where 



-L(d x 9 p ) 2 +K p (d x 4 p ) 2 

9l/2 . 9l/4 AQ 

sin 26 p + cos 40 



H pl =—DV ± ncos29 p + —V^n 2 
na 2a 



(4) 
(5) 



with D = (2 — ^n + Sjj/il + S 2 ) and c 2 w 
[U/(V2~Trt)](a/a) 2 .The. quantity a is a short-length cutoff of 
the order of the lattice constant a. The first term (4) is the 
ID part and takes the same form of the phase Hamiltonian 
for the ID dimerized EHM at quarter- filling, derived in Ref. 
16. 25 Here, the effects of finite lattice dimerization, compared 
with our previous work for 64 = 0, 20 are (i) the appearance of 
the sin 20 p potential and (ii) the modifications of the the cou- 
pling constants (e.g., the factor D in the cos 20 p potential). 
The order parameter n is determined by the self-consistency 
condition 

n = -— D(cos29 p ), (6) 

where (• • • ) denotes the expectation value with respect to H p . 
Equation (6) is obtained from the stationary condition of eq. 
(5) and takes the same form of eq. (14) in Ref. 20 except the 
appearance of the factor D. The equation to determine Tqo is 
given by 



1 = ™LD 2 V ± 



dx 



l/Tco 



dTF CO (x,T)\ (7) 



where F C o{x, r) = (T T cos 29 p (x, t) cos 20 p (O, 0)) is the 
correlation function of the CO order parameter and (• • • ) is 
the expectation value with respect to H p o, i.e., the Hamilto- 
nian of the purely ID system. 

Here we calculate the correlation function based on 
the RG treatment. We also evaluate Fbo( x , t ) = 
(T T sin20 (9 (a;,T)sin20p(0, 0))o, which corresponds to the 
correlation function of the bond-order parameter. Based on 
the conventional RG treatment, these are given by 26 



-Fbo {() = 2 CX P 



- / AC {AK P {1') + 2Gi /4 (0) 
Jo 



d£' (4Kp(C) - 2Gi /4 (0) 



, (8) 



, (9) 




Fig. 2. (Color online) Charge-ordering transition temperature, Tco as a 
function of zV±/t for <5 d = and <5 d = 0.15 for V < V C 1D . The 
inset shows the critical interchain interaction zVf/t as a function of the 
dimerization <5d. 



where £ = ln(r/a) with r = ^Jx 1 + (v p t) 2 . The RG equa- 
tions of K p (£), Gi/±(t), and Gi/ 2 {£) take the same form 
of eqs. (6)-(8) in Ref. 16, with K p (0) = K p , G 1/2 (0) = 
-9i/2/(nv p ), and G 1/4 (0) = g x/ J (2irv p ). In the ID case, 16 
i.e., for V± = 0, the dimer-Mott insulating state is obtained 
for small intersite Coulomb repulsion (V < V^ B when U 
fixed at a large value), where both Giy 2 (^) and Gi^(£) be- 
come relevant under the scaling procedure, while for large in- 
tersite repulsion (V > V^ D ), the COI state is obtained in the 
ground state where Gi/±(£) becomes relevant but Gi/ 2 (£) ir- 
relevant. In the present perturbative RG scheme, the correct 
behavior of the correlation functions is not reproduced for 
large £ due to the relevant behavior of the nonlinear terms 
in the low-energy limit. 26 In order to avoid this problem, we 
simply stop the scaling procedure at the scale corresponding 
to the excitation gap I = For longer distance / > we 
assume that, in the V > V^ 10 case, Fqo CO remams constant 
while -FboCO decays exponentially as a function of r with 
correlation length £ ~ aexp^), and vice versa for -FcoCO 
and F BO W in the V < V C 1D case. 

By inserting the result of eq. (8) into eq. (7) with the above 
procedure, we evaluate Tco- In Fig- 2, we show the transition 
temperature as a function of zVx/t for 5<j = and r5a = 0.15 
in the V < V™ case. For 6a = 0, infinitesimal V± trans- 
forms the system from TLL to COI state when the TLL pa- 
rameter satisfy 1/4 < K p < 1/2, i.e., in the critical region. 20 
The present results show that the dimerization suppresses the 
transition temperature Tco drastically. It is noticeable that a 
finite amount of V± is necessary for the appearance of the 
COI state, for which its critical value V± is shown in the inset 
of Fig. 2. For V± > Vj_, the V± dependence of the transition 
temperature is expressed as 

A-Mott 



T 



CO 



In 



\v±-viJ 



(10) 



where Amou is me Mott gap in the absence of the interchain 
interaction. This formula can be derived by noting that Fco 
decays exponentially. 

In Fig. 3, Tco f° r the V > V^ D case is shown, where, on 
the other hand, the COI state is stabilized even for V± = 0. 
Again, the dimerization suppresses Tco - However, for a fixed 
value of <5d, the critical value Vj_ of the effective interchain 



Fig. 3. (Color online) Charge-ordering transition temperature, Tco as a 
function ofzV±/t for <5 d = and S d = 0. 15 for V > V C 1D . 



interaction V± for the appearance of the CO state is now al- 
ways zero, i.e., infinitesimal interchain interaction leads to fi- 
nite Tco- The V± dependence of the transition temperature 

"1/2 

is now given by Tco fx V, ' as long as Tco *C Aco with 
Aco being the charge gap of the purely ID COI. 

Summarizing the above results, schematic phase diagrams 
on the V±-V^ plane are shown in Fig. 4 where the shaded 
CO regions denote T C o > 0. In the V < V^ D case (a), be- 
tween the anti-phase and in-phase CO states, there appears the 
dimer-Mott insulating phase where CO is absent. On the other 
hand, in the V > case (b), there is a direct transition be- 
tween the two CO states; at the boundary, Tco = 0, i.e., CO 
is realized only at the ground state. 

Next, we focus on the critically fluctuating line, V = V^. 10 . 
The quantum phase transition in the ID case (i.e., V± = 0) is 
known to be the Ising-type, as mentioned. On this line the 
finite-temperature CO phase transition for V± ^ shows 
nontrivial behavior, since the correlation function of the CO 
order parameter Fco t) exhibits a power-law of r -1 / 4 . At 
finite temperature, the analytical form can be obtained as 17 



Fco(x,t) = (T r cos20 p (x,T)cos20 p (O,O)) 



(ttT) 2 



sinh [ttT(x — it)] sinh [irT(x + it)] 



1/8 



(ID 



By inserting this into eq. (7), we find that the V± dependence 

~ 4/7 

of the transition temperature is given by Tco °c Vj_ ■ 

Now we turn to the COI phase below the transition tem- 
perature Tco and derive a useful formula for the relation be- 





Fig. 4. (Color online) Schematic phase diagrams of model (3) on the plane 
of the interchain interactions V± and V[, for the V < V C 1D case (a) 
and for the V > V C 1D case (b). CO and dimer-Mott stand for charge- 
ordered insulating phase with Tco > and dimer-Mott insulating phase, 
respectively. 
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tween splitting in the NMR spectrum 2 ' 3 ' 5 and the CO ampli- 
tude n. In the presence of the Wigner crystal-type CO, an ef- 
fective alternating site potential zV±n(— l) l rn (its wavenum- 
ber is q = it /a = 4fcp) is applied to the system [see, 
eq.(3)] and results in folding of the one-particle dispersion 
relation. Then the Ak-p -component of the density operators 
is expressed in part by the q = density operators around 
the Fermi energy, namely, the slowly-varying component of 
the charge/spin density becomes a part of its 4/c F compo- 
nent. This is directly related to the splitting in the NMR spec- 
trum since the response under uniform magnetic field at each 
site becomes disproportionated with a period of two sites 27 
reflecting the Knight shift. In the phase representation, the 
staggered component of the spin density m^,, (x) is propor- 
tional to the slowly varying component of the spin density 
m (x) = d x 9 a (x)/ir as 



m 4 fc F (x) = -(-!)' 



zV ± n 



d x 9 a (x) 



(12) 



From this relation, the Knight shift at the ith site can be ex- 
pressed as 



Si(X)cxXa(T) 



l-(-l)' 



zV±n 



^2t^l + 5l) + {zV^nf 



(13) 



where \a (T) is the magnetic susceptibility. This verifies that 
the average of the Knight shift is proportional to the mag- 
netic susceptibility as observed experimentally. 2 Furthermore 
we find that the "relative shift" (Si - S )/(Si + S ) can be 
expressed as 



Si — So 
Si + So 



zV±n 



2t 2 (l 



(14) 



{zV^nf 



We note that the intrachain interactions do not explicitly ap- 
pear in this expression, but they are embedded in the mag- 
nitude of the order parameter n. Experimentally, the order 
parameter of the CO phase has been sometimes estimated 
so far by simply assuming that the Knight shift is propor- 
tional to an amount of the charge located on the site, i.e., 
(Si - So) /(Si + S ) — > 2n. 2 However, the value of n ob- 
tained in such a way can have been misestimated. 

We compare the present results with recent experimental 
observations by dielectric permittivity 28 and ESR 29 that Tco 
for some (TMTTF-<ii2)2 A" salts composed of the perdeuterio- 
TMTTF are increased. The X-ray analysis for the TMTTF- 
d\2 and TMTTF-/112 compounds has revealed that the deuter- 
ation leads to suppression of the lattice dimerization, 29 which 
is consistent with the present results. The rather large rise in 
Tco °f about 10 percent compared to the small change in 
dimerization suggests that the systems are in the critical re- 
gion, V < V C 1D . 

The anomalous behavior of the dielectric constant, pro- 
portional to \T — T co | _1 , observed in the TMTTF com- 
pounds, 430 has been theoretically discussed based on the 
dimerized EHM in the purely ID case. 4, 18 It is suggested that 
at T = the dielectric constant shows a divergent behavior 
at the boundary between the Mott insulator and the COI, i.e., 



V = V^ 10 . If the interchain coupling is taken into account 
as in our work here, this quantum phase transition would be 
changed into a finite-temperature transition and we expect that 
the dielectric constant should show a similar anomalous be- 
havior at T"co- The investigation of such a possibility in the 
present scheme remains for a future study. 

In conclusion, we have discussed the competition between 
the charge-ordered insulating state and the dimer-Mott insula- 
tor at finite temperatures in a quasi-one-dimensional extended 
Hubbard model at quarter filling. Frustration in the interchain 
Coulomb repulsion reduces the effective interchain interac- 
tion and therefore decreases the charge-ordering transition 
temperature Tco- The lattice dimerization also reduces Tco, 
while a finite strength of the effective interchain interaction 
is necessary to transform the system from the dimer-Mott in- 
sulating state into the charge-ordered insulating state. Based 
on the present analysis, we have derived a useful formula for 
estimating the magnitude of the order parameter of the CO 
phase experimentally. 
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